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ABSTRACT 

We report one of several homologous non-radial eruptions from NO A A active region (AR) 11158 
that are strongly modulated by the local magnetic field as observed with the Solar Dynamic Observa- 
tory (SDO). A small bipole emerged in the sunspot complex and subsequently created a quadrupolar 
flux system. Non-linear force-free field (NLFFF) extrapolation from vector magnetograms reveals its 
energetic nature: the fast-shearing bipole accumulated ~2xl0 31 erg free energy (10% of AR total) 
over just one day despite its relatively small magnetic flux (5% of AR total). During the eruption, 
the ejected plasma followed a highly inclined trajectory, over 60° with respect to the radial direction, 
forming a jet-like, inverted-Y shaped structure in its wake. Field extrapolation suggests complicated 
magnetic connectivity with a coronal null point, which is favorable of reconnection between different 
flux components in the quadrupolar system. Indeed, multiple pairs of flare ribbons brightened simul- 
taneously, and coronal reconnection signatures appeared near the inferred null. Part of the magnetic 
setting resembles that of a blowout- type jet; the observed inverted-Y structure likely outlines the 
open field lines along the separatrix surface. Owing to the asymmetrical photospheric flux distribu- 
tion, the confining magnetic pressure decreases much faster horizontally than upward. This special 
field geometry likely guided the non-radial eruption during its initial stage. 

Subject headings: Sun: activity — Sun: corona — Sun: surface magnetism — Sun: magnetic topology 



1. INTRODUCTION 

Solar eruptive events derive their energy from the non- 
potential coronal magnetic field (Forbes 2000; Hudson 
2011). Reconnection takes place locally where the field 
gradient is large, but can alter the larger-scale field topol- 
ogy rapidly. The dissipated energy from the relaxing field 
accelerates particles, produces radiation, and heats and 
ejects plasma into the interplanetary space as a coronal 
mass ejection (CME). 

Prior to eruption, energy builds up in the corona 
through flux emergence and displacement, which may 
take up to a couple of days (Schrijver 2009). The 
slow evolution can be approximated by a series of 
quasi-stationary, force-free states in the low plasma- 
f3 coronal environment. This allows the estima- 
tion of AR energetics in non-flaring states, thanks 
to recent advances in photospheric field measurement 
and field extrapolation algorithms (Regnier & Canfield 
2006; Thalmann & Wiegelmann 2008; Jing et al. 2009; 
Sun et al. 2012). 

Besides the gross energy budget, the detailed magnetic 
configuration also proves important to the initiation, ge- 
ometry, and scale of eruptions. In the case of a coro- 
nal jet, the direction of the ambient field (horizontal or 
oblique) directly determines the direction of the jet and 
its distinct emission features (two-sided or "anemone" 
type, Shibata et al. 1997). Observation and modeling 
demonstrate that the overlying field provides a critical 
constraint on CME's speed and trajectory (Liu 2007; 
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Gopalswamy et al. 2009; Wang et al. 2011). 

Theoretical studies have extensively explored the role 
of topological features in reconnection (Demoulin et al. 
1996; Priest & Forbes 2000; Longcope 2005). Their 
applications to solar events usually involved the 
results of potential or linear force-free field ex- 
trapolation (Aulanier et al. 2000; Fletcher et al. 2001; 
Mandrini et al. 2006), or magnetohydrodynamic (MHD) 
simulations that qualitatively reproduce the observed 
phenomena (Moreno-Insert is et al. 2008; Pariat et al. 
2009; Masson et al. 2009; Torok et al. 2011). 

Here we report one of several similar non-radial erup- 
tions that are strongly modulated by the local mag- 
netic field as observed with the Solar Dynamic Observa- 
tory (SDO). Using vector magnetograms from the Helio- 
seismic and Magnetic Imager (HMI) (Schou et al. 2012; 
Hoeksema et al. 2012) aboard SDO and a non- linear 
force-free field (NLFFF) extrapolation, we monitor the 
AR evolution and explain the magnetic topology that 
leads to the curious features during the eruption. The At- 
mospheric Imaging Assembly (AIA; Lemen et al. 2012) 
and other observatories recorded these features and pro- 
vide guidance for our interpretation. 

In Section 2 we briefly describe the data and the ex- 
trapolation algorithm. We first present observations of 
the eruption in Section 3, and then come back in Sec- 
tion 4 to explain the magnetic field and energy evolution 
leading to the event. In Section 5, we interpret this cu- 
rious event based on the magnetic field topology. We 
discuss the interpretation in Section 6 and summarize in 
Section 7. 
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Fig. 1. — Full-disk, unsharp masked AIA 171 A image at 17:28:15 
UT on February 14, 2011 showing the non-radial eruption. Inset 
shows the enhanced image of the ejecta. The two flux-rope-like 
structures with a shared eastern footpoint are marked as FR1 
and FR2. Animation of a 20-hr interval shows at least five 
similar eruptions. (An animation of this figure is available at 
http: //sun. Stanford. edu/~xudong/Article/Cusp/homolog.mp4.) 



Sunspot complex AR 11158 produced the first X-class 
flare of cycle 24 near its center on 2012 February 15 
(Schrijver et al. 2011). Before and after that flare, there 
were a series of smaller eruptions from its northeastern 
periphery, our region of interest (ROI), where a small 
new bipole emerged. Five of them assumed very similar 
structures and were accompanied by C or M-class flares 
within a 20-hr interval (06:58, 12:47, 17:26, and 19:30 UT 
on February 14, and 00:38 UT on February 15; see the 
animation of Figure 1 and Figure 4(d)). In all cases the 
ejecta followed a similar, non-radial trajectory towards 
the northeast. 

We focus here on the event around 17:26 UT on Febru- 
ary 14 associated with an M-2.2 class flare. The eruption 
site was near central meridian (W04S20). For context, 
we study the AR field evolution during a 36-hr interval 
leading to and shortly afterward the event, from Febru- 
ary 13 12:00 UT to February 15 00:00 UT. 

The HMI vector magnetograms provide photospheric 
field measurement at 6173 A with 0.5" pixels and 12- 
minute cadence. Stokes parameters are first derived from 
filtergrams averaged over a 12-minute interval and then 
inverted through a Milne-Eddington based algorithm, 
the Very Fast Inversion of the Stokes Vector (VFISV; 
Borrero et al. 2011). The 180° azimuthal ambiguity in 
the transverse field is removed using an improved ver- 
sion of the "minimum energy" algorithm (Metcalf 1994; 
Leka et al. 2009). Here, the selected 36-hr dataset in- 
cludes 181 snapshots of a ^OO^xSOO" region. For data 
reduction procedures, we refer to Hoeksema et al. (2012) 
and references therein. 

We use an optimization-based NLFFF extrapolation 
algorithm (Wiegelmann 2004) and HMI data as the lower 
boundary to compute the coronal field. The side and 
upper boundaries are determined from a potential field 



extrapolation (PF) using the Green's function method 
(Sakurai 1989). The computation domain assumes pla- 
nar geometry, uses a Cartesian grid (300x300x256) 
and a 720 km (~1") resolution. Before extrapola- 
tion, we apply to the data a pre-processing procedure 
(Wiegelmann et al. 2006) that iteratively reduces the net 
torque and Lorentz force so the boundary is more con- 
sistent with the force-free assumption. The magnetic 
free energy is simply the energy difference between the 
NLFFF and PF. Our previous study on the same region 
(Sun et al. 2012) used identical procedures, where we de- 
scribed and evaluated the algorithm in detail. 

3. THE NON-RADIAL ERUPTION 

Observed in the AIA extreme- ultraviolet (EUV) bands, 
a small AR filament situated above the polarity inversion 
line (PIL) of a newly emerged bipole started its slow rise 
around 17 UT (see online animation of Figure 2). The M- 
class flare peaked at 17:26 UT in soft X-ray (SXR) flux, 
when the filament rapidly erupted towards the northeast. 
The ejecta appeared to consist of two rope-like features 
(FR1 and FR2 in Figure 1) with a shared eastern foot- 
point. By inspecting AIA image sequences in various 
bands and HMI magnetograms, we think that they orig- 
inated from the same filament structure. 

The STEREO-A spacecraft was then near quadrature 
with the Sun-Earth line (87° ahead). Its SECCHI EUVI 
instrument (Howard et al. 2008) caught a glimpse of the 
ejecta in the 195 A channel (Figures 2(a)), where the 
erupted filament appeared to follow a straight trajectory 
viewed from west. Using the simultaneous image from 
the AIA 193 A channel (Figures 2(b)), we are able to 
estimate its three-dimensional (3D) geometry. 

Figure 2(c) illustrates the triangulation procedure. We 
manually select the eruption site O and the frontmost 
point P of the inner flux rope FR2 (as projected on the 
plane of sky) in the AIA image. We select the corre- 
sponding points O' and Q in the EUVI image, such that 
1) O and O' have the same Carrington coordinate; 2) the 
ejecta's N-S extent in two images satisfies |OR| = JO'R'I, 
where OR and O'R' represent the projection of line seg- 
ments OP and OQ in the N-S direction in SDO's plane- 
of-sky, respectively. 

Assuming the ejecta follows a straight trajectory, we 
can solve for its inclination 5 and azimuth a with respect 
to the line-of-sight (LOS). We find that (5=43°, a=34°. 
By repeating the point selection process we estimate the 
uncertainty to be ~3° under the current scheme. The 
trajectory is highly inclined, about 66° with respect to 
the local radial direction. 

A bright, inverted- Y shaped structure formed in the 
wake of the eruption. It consisted of a thin spire on 
top of a cusp-shaped loop (Figure 2(d)); both lasted 
over 1 hr. The cusp appeared almost two-dimensional 
and had both "legs" rooted in negative polarity flux 
(see Section 5 and Figure 5(b)). There were propagat- 
ing brightness disturbances along the cusp legs and the 
spire (Thompson et al. 2011, see the animation of Fig- 
ure 2), which have been interpreted as episodic plasma 
flows (see the coronal seismic and Doppler analyses in 
Su et al. 2012; Tian et al. 2012). These observed features 
outline a magnetic arrangement that resembles a coronal 
jet (e.g. Shibata et al. 1997). Nevertheless, the struc- 
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Fig. 2. — Geometry of the non-radial eruption. O and (V mark the eruption site, (a) SECCHI EUVI 195 A image from STEREO-A, 
about 87° ahead of SDO. Due to the tilt of the solar rotational axis, the SDO and STEREO north are offset by 6.8°. (b) AIA 193 A image 
of the same ejecta, taken 4 s later than (a). The projected N-S length of the ejecta (|OR|) is identical to that in (a) (IO'R/1), where OR 
and / R / represent the projection of line segments OP and OQ in the N-S direction in SDO's plane-of-sky, respectively. The scales of (a) 
and (b) are different in order to better show the features of interest, (c) Schematic diagram explaining the determination of the ejecta's 
geometry. SDO's west, north, and LOS directions are taken as x, y, and z axis. The pink arrow represents the ejecta, its projected shape 
viewed from EUVI and AIA are shown as pink dashed lines on green and brown planes. The local radial vector is about W13S04 to LOS. 
The inclination 5 is about 43°; azimuth a about 34°. See Section 3 for details, (d) AIA 171 A image of the post-eruption AR; Y marks 
the top of the cusp and the base of the jet. The boxed region is used to construct panel (e). Purple/pink contours are for HMI LOS field 
at ±200 G. (e) Space-time diagram showing the speed of ejecta and jet. Three dashed lines (starting near 17:25, 17:36, and 17:49 UT) 
indicate a projected speed of 500, 330 and 280 km s _1 , respectively. Panels (a), (b), and (d) are displayed in a square-root scale. (An 
animation of this figure is available at http://sun.stanford.edu/~xudong/Article/Cusp/ejecta.mp4.) 



ture appeared only after the eruption. Various observed 
features appear to require alternative explanations other 
than the standard jetting model or its variations (see a 
brief discussion in Section 6.3). 

By placing a cut along the thin spire in the AIA 171 
A image sequence, we construct a space-time diagram to 
illustrate the relevant speeds in this event (Figure 2(d) 
and (e)). The projected speed of the ejecta is about 
500 km s _1 ; the brightness disturbance is around 300 
km s _1 . Considering the inclined trajectory, we estimate 
the real speed about 30% higher, i.e. 650 and 390 km 
s _1 , respectively. 

4. THE EMERGING BIPOLE AS ENERGY SOURCE 

We study the underlying photospheric field that led 
to this eruption. Figure 3(a) shows a snapshot of the 
radial field taken 25 minutes before the event as derived 
from the vector magnetogram. The AR mainly consists 
of two interacting bipoles. A large amount of magnetic 
free energy was stored near the major PIL between the 
shearing sunspots at center of the field of view (FOV), 



where the X-class flare took place (Sun et al. 2012). 

The eruption studied here is related to a newly 
emerged, smaller bipole (boxed region in Figure 3(a)). 
The bipole appeared on February 13 in the northeast- 
ern part of the AR. Starting from 12 UT on February 
14, the positive component advanced rapidly westward 
with strong rotational motion and shearing with respect 
to its negative counterpart, leaving behind a fragmented 
stripe of flux mimicking a long-tailed tadpole (see the 
online animation of Figure 3). 

The new bipole had strong horizontal photospheric 
field that that lay parallel to the PIL (Figure 3(b)). 
NLFFF extrapolation suggests a highly twisted core field 
and strong radial current (Figure 3(c)), which correspond 
to the observed AR filament that eventually erupted. 

We summarize in Figure 4 the bipole's temporal evo- 
lution. By evaluating the area within the ROI (boxed 
region in Figure 3(a)), we estimate its unsigned flux to 
be only about 5% of the AR's total around the eruption 
time. However, the surface unsigned radial current with 
in the ROI accounts for 12% of the AR's total, much 
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Fig. 3. — Snapshot of magnetic field of AR 11158, 25 minutes before the eruption, (a) Radial magnetic field (B r ) map as derived from the 
vector magnetogram. The contours are for ±200 G. PI, Nl, P2, and N2 mark four components of the quadrupolar flux system. The yellow 
box indicates the FOV for (b) and (c) and is identical to that in (d). (b) Photospheric vector magnetic field map. Gray-scale background 
shows B r . The blue/red arrows indicate the horizontal component (B^) with positive/negative radial counterpart, where field strength 
B > 200 G. Their lengths correspond to the magnitude (B^)] their directions show the azimuth, (c) Selective extrapolated field lines plotted 
on B z map. The color shows the amount of radial current at the field line footpoint. (d) Map of current density (| J\) integrated over the 
lowest 10 Mm in extrapolated field. The light/dark gray contours are for B r = ±200 G. All data are deprojected and remapped using the 
Lambert equal area projection. (An animation of this figure is available at http://sun.stanford.edu/~xudong/Article/Cusp/field.mp4.) 
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Fig. 4. — Evolution of the entire AR and the region of interest 
(ROI). (a) Unsigned magnetic flux, (b) Unsigned photospheric 
electric current, (c) Magnetic free energy evaluated in the volume 
above the ROI. (d) GOES soft X-ray flux. The homologous flares 
are marked with "M" and "C" according to their GOES X-ray 
class. In panels (a)-(c), solid line represents the whole AR. Dotted 
line is for the bipole (boxed region in Figure 3(a)); bipole values are 
multiplied by 5 for clarity. The gray vertical band indicates a 2.6- 
hr period that brackets the eruption, whose means are compared 
with the AR means. For (a) (b) only pixels with B > 200 G 
are included. Errors in (a) (b) derived from spectropolarimetry 
inversion are small. Errors (3<r) in panel (c) are evaluated using 
a pseudo Monte-Carlo method (Sun et al. 2012) which show the 
effect of spectropolarimetric noise. 



higher than the corresponding flux fraction. We inte- 
grate the free energy in the volume above the ROI and 
find it to be over 10% of that in the whole volume. For 
the ROI, the ratio between the NLFFF energy and the 
PF energy is about 1.60. This indicates the bipole is very 
non-potential and energetic. There is a strong concentra- 
tion of current near the PIL in the lower corona, similar 
to the major PIL near the center of the AR (Figure 3(d)). 

Unfortunately, we do not find a clear, step- wise change 
in free energy during the flare that can be used as a proxy 
of the energy budget (Figure 4(c)). Our earlier work on 
the ensuing X-class flare (Sun et al. 2012) suggests the 
energy budget tends to be underestimated by the extrap- 
olation method. This is partly because the flaring field is 
dynamic and likely not force- free (e.g. Gary 2001); thus 
it cannot be reliably described by the NLFFF model. 
Limited resolution and uncertainties in the field mea- 
surement and modeling may also be a factor. The free 
energy for the ROI gradually decreased after 20 UT when 
the positive flux fragmented and the current decreased. 

The emergence of the bipole led to a local enhancement 
of free energy with a series of ensuing eruptions from 
this relatively small region. Its very existence changed 
the original magnetic configuration and converted it into 
an asymmetrical (the new bipole is relatively small) 
quadrupolar flux system. The change of the photospheric 
flux distribution altered the coronal magnetic connectiv- 
ity in a fundamental way, and may have contributed to 
the destabilization of the system. For clarity, we label 
the four quadrupolar components PI, Nl (including the 
old sunspot and the negative part of the new bipole), P2, 
and N2 (Figure 3(a)). 
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Fig. 5. — Magnetic topology based on NLFFF extrapolation for the pre-emption state, (a) SDO view of four sets of loops connecting the 
four quadrupolar flux components pairwise, as well as twisted field lines below representing the AR filament. The cross section is identical 
to that in (d). Inset shows the corresponding AIA 94 A image, which is the same as Figure 6(a). The inferred coronal null point, marked by 
"X", appears slightly above the observed loops, (b) Magnetic null point, spine field line, and open field lines that outline the separatrix (fan) 
surface, (c) Side view of the region (from east), (d) Side view with z-axis (radial direction) stretched by 2. Magnetic pressure is imaged on 
a vertical cross section to illustrate its anisotropy. The cross section is roughly aligned with the direction of eruption, and is in front of the 
null from this viewing angle. (An animation of this figure is available at http://sun.stanford.edu/~xudong/Article/Cusp/topo.mp4.) 



5. INTERPRETATION BASED ON THE 
MAGNETIC FIELD TOPOLOGY 

5.1. A Coronal Null and the Inclined Trajectory 

What is the coronal magnetic field topology that led 
to the highly non-radial eruption? Field lines computed 
from the pre-flare NLFFF solution (16:59 UT) reveal con- 
nectivity between each pair of the opposite polarity flux 
(Pl/Nl, P2/N1, P2/N2, and P1/N2) in this quadrupo- 
lar system (Figure 5(a)). Such connectivity is apparent 
in the AIA observations. 

One striking feature, however, is the large gradient 
in field line mapping. For example, loops connecting 
P2/N1 (cyan) and P2/N2 (orange) are at first parallel, 
but diverge drastically near their apexes, becoming al- 
most antiparallel with each other. These modeled field 
lines closely resemble the observed loops (inset of Fig- 
ure 5(a) or Figure 6(a)). The cusp-like P2/N2 and the 
diverging field lines strongly suggest the existence of a 
coronal null point, where field strength becomes zero. 3 

Using a trilinear method (Haynes & Parnell 2007), we 
indeed find a null point situated at ~9 Mm height (Fig- 
ure 5(b)) right above the modeled loop apexes (see Ap- 
pendix A). From that null, closed loops "turn away" with 
a sharp angle. Seen from side (Figure 5(c)), these loops 
are low- lying; they incline towards the northeast, the di- 
rection of the eruption. This configuration persisted over 

3 See TRACE observation of AR 9147/9149 
(http : //trace . lmsal . com/P0D/TRACEpodarchive4 . html ). 



the next few hours (see Section 6.1). 

This inclined geometry is perhaps a natural conse- 
quence of the asymmetrical photospheric flux distribu- 
tion. We infer that this field configuration may have 
facilitated the non-radial eruption in the following ways. 
First, reconnection may take place near the null point, 
removing the overlying flux above Pl/Nl and preferen- 
tially reducing the confinement from the northeast direc- 
tion. Second, the ambient, confining magnetic pressure 

(Pb—§^) is anisotropic: it drops off much faster horizon- 
tally than it does in the radial direction (Figure 5(d)). 
When the anisotropy is strong enough, it can guide the 
ejecta towards a direction with large negative pressure 
gradient by deflecting its trajectory. It effectively cre- 
ates a non-radial "channel" for the plasma to escape. 

5.2. The Inverted- Y Structure 

We further analyze the magnetic topology of the pre- 
emption state for insight on the observed inverted-Y 
shaped structure. By analyzing the Jacobian field matrix 
(Mij = dBi/dxj) at the inferred null point, we are able 
to find the spine and the fan, which are special field lines 
that define the magnetic configuration near the singular- 
ity (e.g. Parnell et al. 1996). Regular field lines passing 
by the immediate vicinity of the null point generally out- 
line the separatrix (fan) surface (Figure 5(b)(c)). In this 
case they separate the closed flux inside and the open 
flux outside. We describe the analysis method in Ap- 
pendix A. 
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Fig. 6. — Various flare emission observations. The time of each image is marked in (g). (a) Composite of negative AIA 94 A image (~6 
MK, 2-minute average) and HMI LOS magnetogram showing the possible coronal null point. The AR filament is manually outlined with 
dotted line. The null is at (77" , -209" ) on the plane of sky. See also Figure 5(a). (b) Hinodej 'SOT Ca II H band showing flare ribbons 
(Rpi, Rp2, Rni 5 R-N2) and the underlying photosphere. Half-ring-like secondary ribbon (Rno) is a ls° visible, (c) Unsharp masked, negative 
SOT Ha image showing flare ribbons and magnetic connectivity. An additional, weak brightening is marked as Rpo- The erupting plasma 
is visible in the foreground from about (60", -250") upward; the tip of the cusp is near (45", -160"). (d) RHESSI 12-25 keV HXR image 
as contours on negative AIA 131 A image (~10 MK), showing a footpoint source (SI) and coronal source (S2). The two ribbons Rpi 
and Rni are spatially unresolved in HXR. Contours are drawn at 20%, 50%, and 90% of maximum, (e) Same as (d), for 4 minutes later. 
Contours are for 10% and25% of the maximum of (d). (f) Hinode /XRT SXR image, showing the cusp-like structure, (g) GOES SXR flux, 
Fermi /GBM 12-15 keV, and 25-50 keV HXR flux. RHESSI coverage of the flare started from 17:27:44 UT, and is not shown here. Panels 
(b)-(e) are displayed using square-root scale. All images are tracked with solar rotation and co-aligned to an accuracy better than 0.6". 



Owing to the local excess of negative flux, open field 
lines from Nl and N2 flow along the separatrix and con- 
verge around the outer spine. These field lines naturally 
form an inverted- Y structure (Figure 5(b)). Its morphol- 
ogy resembles the observed loops, although less inclined 
towards the northeast. Their detailed geometry took 
shape during the the dynamic eruption, which the static 
extrapolation is unable to model. 

5.3. Observational Evidences 

Because field line mapping diverges and links the whole 
quadrupolar system, we expect electrons accelerated dur- 
ing the flare near the null point to precipitate along differ- 
ent loop paths, resulting in multiple pairs of flare ribbons 
brighting simultaneously (Shibata et al. 1995). Taken by 
the Solar Optical Telescope (SOT; Tsuneta et al. 2008) 
on the Hinode satellite, Ca II H band images (Fig- 
ure 6(b)) indeed show such phenomena. The typical 
double ribbons (Rpi/Rni) are related to the erupting 



filament, whereas Rp2 and Rn2 are likely related to the 
reconnecting P2/N2 loop. Ka images (Figure 6(c)) pro- 
vide additional information on the magnetic connectivity 
between Rpi/Rn2 and Rp2/Rni- Remarkably, the rib- 
bon Rp 2 appears to be co-spatial with the inferred spine 
field line footpoint (Figure 5(b)), which moved with time 
as seen in the Ca II H and Ha image sequences. 

The Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI] Lin et al. 2002) missed the impulsive phase 
but captured what appeared to be a coronal hard X- 
ray (HXR) source (S2 in Figure 6(d)) in the flare's early 
decaying phase. The source's proximity to the inferred 
coronal null gives strong support to our interpretation. 
From the loop top, energetic electrons followed very in- 
clined paths towards the footpoints in Pl/Nl, which 
created the footpoint source (SI) corresponding to the 
Rpi/Rni ribbons. This coronal source lasted well into 
the decaying phase (Figure 6(e)). 
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6. DISCUSSION 

6.1. On the Coronal Field Topology 

How common is the magnetic topology determined 
here? A previous study focused on the quadrupolar 
configuration of AR 10486 during the 2003 X-17 flare 
(Mandrini et al. 2006). The major eruption was found 
to involve reconnection at the quasi-separatrix layers 
(QSL; Demoulin et al. 1996), while a smaller brighten- 
ing was associated with a similar coronal null point de- 
termined using a linear force-free extrapolation. In an- 
other quadrupolar region AR 11183, similar cusp and jet 
structures existed at a much larger scale (Filippov et al. 
2012). The white-light jet extended over multiple solar 
radii. 

We analyze the entire 36-hr series, searching for con- 
sistency in time. The coronal null at 9 Mm appeared 
in a few frames early on February 14, distinct from all 
other candidates which were mostly below 4 Mm in weak 
field regions. Starting from 15:35 UT, it appeared at a 
nearly constant location (within 3 Mm of the first de- 
tected null) in over half the frames afterwards (22/42, 
until February 15 00:00 UT), while the near-surface nulls 
rarely repeated in two consecutive time steps. We have 
applied a different null-searching method based on the 
Poincare index theorem (Greene 1992) and found similar 
results (23/42, 20 identical to the trilinear method, with 
3 additional and 2 missed detections). The repeated de- 
tection of null points and the observed homologous erup- 
tions (Figure 4(d)) suggest the aforementioned topology 
is characteristic for this quadrupolar system. 

We compute at 1-hr cadence the "squashing factor" 
Q that describes the field mapping gradient (Titov et al. 
2002) by tracing individual field lines and measuring the 
differences between the two footpoint locations. High-Q 
isosurface corresponds to QSLs. By inspecting the con- 
tour of Q at different heights, we find that multiple QSL's 
tend to converge and intersect at about 9 Mm. Near 
the intersection, the field strength is weak, and the field 
line mapping gradient is invariably large, with or without 
null point. This illustrates the robustness of our inter- 
pretation despite the uncertainties in the extrapolation 
algorithm (e.g. DeRosa et al. 2009) and the field mea- 
surement. (The uncertainties nevertheless can indeed af- 
fect the detailed fan-spine configuration, as discussed in 
Appendix A.) 

We note that our PF extrapolation, with radial field 
as boundary condition and the Green's function method, 
does not detect any nulls above 5 Mm. Instead, we find 
a low- lying null at about 4 Mm in 13 frames, southwest 
to the NLFFF solution. The field configuration is less re- 
alistic, presumably because the current-free assumption 
does not agree with observation. 

6.2. On the Flare Emissions 

Owing to the LOS projection, the altitude of an on-disk 
HXR source cannot be unambiguously determined. We 
think S2 is a coronal source mainly because it appeared 
near the apex of cusp-shaped loops (P2/N2) which is typ- 
ical for reconnecting field lines (e.g. Tsuneta 1996). In 
addition, its strong HXR emission (peak at ~60% of the 
maximum) does not correspond to any bright flare rib- 
bon. The closest chromospheric emission enhancement is 
a small patch (Rpo in Figure 6(c)) within a fragmented 




Fig. 7. — Schematic illustration of the magnetic configuration 
and dynamics that may have led to the eruption. The structure 
resembles that of a blowout jet. The arcade (blue field lines above 
Pl/Nl) from the newly emerged bipole expands, reconnects with 
the pre-existing field (blue field lines from N2), becomes open (yel- 
low field lines from Nl), and the low- lying sheared/twisted core 
field (pink field lines between Pl/Nl) subsequently erupts. A pos- 
sible initial reconnection site near is marked by the star; possible 
motions of the loops are denoted by thick arrows. Pre- and post- 
reconnection field lines are colored blue and yellow, respectively. 
The directions of the observed, post- eruption flow (Figure 2 and 
animation, see also Thompson et al. (2011); Su et al. (2012)) are 
denoted by thin arrows. The inset shows the SXR difference image 
between 17:22:32 and 17:19:56 UT from Hinode XRT Ti Poly fil- 
ter (FOV 72"x60"). The brightening P1/N2 loop is marked by a 
yellow circle; the brightening filament is visible in the foreground. 



positive flux about 5" to the east and south, whose in- 
tensity is much weaker than the Rpi/Rni ribbons. This 
argues against the footpoint source interpretation. 

We notice a dimmer, half-ring-like ribbon (Rno) far- 
ther north in the weak field area (Figure 6(b)); both 
Ha (Figure 6(c)) and EUV images (animation of Fig- 
ure 2) show its connection to PI. This structure is 
related to flux emerging into an encircling unipolar re- 
gion ("anemone" AR; Shibata et al. 1994). Because the 
brightening Rno region possesses flux only a few percent 
of PI (c.f. Reardon et al. 2011), we consider this struc- 
ture secondary. It does not affect our conclusions on the 
AR topology. 

Because no HXR source was detected at the P2/N2 
footpoints and the Rp2/Rn2 ribbons were fainter than 
Rpi/Rni, we think the electrons primarily precipitated 
along the shorter Pl/Nl loop during the flare. On the 
other hand, the P2/N2 loop produced much stronger 
SXR and EUV emission during the flare's late de- 
caying phase. Almost 30 minutes later, SXR images 
(Figure 6(f)) from the Hinode X-Ray Telescope (XRT; 
Golub et al. 2007) still showed a bright cusp structure 
above P2/N2. 

6.3. On the Eruption Mechanism 

When a bipole emerges, one leg of the new loop may 
reconnect with the oppositely directed, pre-existing open 
field. The released magnetic energy heats the plasma and 
produces field-collimated outward flows, known as the 
"standard" jet phenomenon (Shibata et al. 1997). When 
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the emerging field is sheared or twisted, its core may sub- 
sequently erupt. Events in this sub-class have recently 
been described as "blowout" jets (Moore et al. 2010). 

Can this event be explained by the jet models? We 
find the inferred magnetic structure here resembles the 
blowout type. Illustrated in Figure 7, the newly emerged 
bipole (PI and the north part of Nl) hosts a twisted 
core field. We speculate that the increasing flux leads to 
the expansion of the arcade loops above, which reconnect 
with the open, negative-polarity field from N2. This pro- 
cess opens up the arcade loops and acts to promote the 
eventual eruption of the core field below. The jet model 
predicts the brightening of the reconnected P1/N2 loop, 
which is indeed observed in the SXR images (inset of 
Figure 7). However, in contrast to the expected jet be- 
havior, no outward flows are observed during this stage. 
The jet-like, inverted-Y structure appeared only after the 
core field eruption and the accompanying M-class flare. 

Propagating brightness disturbances in the post- 
eruption inverted-Y structure have been interpreted as 
pulsed plasma flow (Su et al. 2012; Tian et al. 2012). 
The upflow from the left leg diverges and flows in oppo- 
site directions, upward in the thin spire and downward in 
the right leg (Figure 7 and Animation of Figure 2). The 
flow is most pronounced in cooler EUV wavebands (e.g. 
171 A, ~0.6 MK) and is absent in SXR images. In the 
standard jet model, these collimated flows are produced 
and heated by reconnection. The relatively low temper- 
ature observed here suggests a low-altitude reconnection 
site with cooler plasma supply (c.f. Su et al. 2012), rather 
than the one near the base of the spire higher in the 
corona. The detailed dynamics of this event require fur- 
ther investigation which is out of the scope of this work. 

7. SUMMARY 
We summarize our findings as follows. 

— Bipole emergence and shearing in a pre-existing 
sunspot complex introduced a large amount of free 
energy, despite its small flux. The new flux pow- 
ered a series of homologous, non-radial eruptions. 

— One typical eruption had an inclined trajectory 
about 66° with respect to the radial direction. 



An inverted-Y structure consisted of cusp and jet 
formed in the wake of the eruption. 

— The bipole emergence created an asymmetrical 
quadrupolar flux system. Field extrapolation sug- 
gests that the consequent, inclined overlying loops 
and the anisotropic magnetic pressure are respon- 
sible for the non-radial eruption. 

— Extrapolation suggests a coronal null point at 
about 9 Mm, slightly below the apexes of the cusp- 
like loops. Its location is favorable for reconnection 
between different flux components in the quadrupo- 
lar system. The observed inverted-Y structure is 
likely related to the open negative field lines in part 
outlining the separatrix surface. 

— Multiple flare ribbons brightened simultaneously 
during the accompanying flare. A coronal HXR 
source appeared near the inferred null point. These 
observations support our interpretation. 

— The inferred magnetic structure resembles that of a 
blowout-type jet. Some observed features fit in the 
jet model, while others remain difficult to explain. 

The event studied here demonstrates the importance of 
detailed magnetic field topology during solar eruptions. 
Flux emergence in suitable environment can lead to fun- 
damental changes in the coronal field geometry, which 
then place strong constraints on the plasma dynamics. 

We thank B. J. Thompson for bringing this event to 
our attention and the anonymous referee for the helpful 
comments. We are grateful to T. Wiegelmann for provid- 
ing the NLFFF extrapolation code. We benefited from 
discussions with M. Derosa, W. Liu, C.-L. Shen, and L. 
Tarr. The SDO data are courtesy of NASA and the 
HMI and AIA science teams. We acknowledge the use of 
STEREO /SECCm EUVI, Hinode/ SOT, XRT, RHESSI, 
GOES and Fermi /GBM data. Figures 5 and 7 are pro- 
duced by VAPOR (www.vapor.ucar.edu). 

Facilities: SDO, Hinode, STEREO, RHESSI, GOES, 
Fermi. 



APPENDIX 

A. METHOD FOR FINDING THE MAGNETIC TOPOLOGICAL SKELETON 

At a magnetic null point, field strength becomes 0, and singularity arises. We follow the null-searching method 
described in Haynes & Parnell (2007). Assuming the field is trilinear within each volume element, the 3D field vector 
B = (Bi, £?2, Bs) T and its derivatives (dBi/dxj, i,j = 1, 2, 3) within each cell are completely determined by the values 
on its eight vertices. To search for possible null point, we first scan over each cell in the domain: if any IV s have the 
same sign on all the vertices, the cell cannot host a null point and will be ignored. For each remaining cell, we use a 
Newton-Raphson scheme to iteratively solve for x = (xi,X2^xs) T that satisfies £^(x) = 0: 

£P) V n )> ( ai ) 

where n and n + 1 denote two consecutive iteration steps, and the repeated index j means summing of all j's. For the 
16:59 UT frame, we find a null point at x = (89.2425, 173.8625, 12. 7667) T in the (300x300x256) domain. At a 720 
km resolution, its height is about 9.2 Mm. The field strength \B\ is about 10 -5 G. 
The rest of the method description is adapted from Parnell et al. (1996) and Haynes & Parnell (2010). To first order, 
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the magnetic field near a null point located at x' is approximated by 

Bi = Mij{xj — Xj), (A2) 

where the matrix M^ = dBi/dxj is the Jacobian matrix, and is evaluated in this case as 

/ dB 1 /dx 1 dB 1 /dx 2 dB 1 /dx 3 \ ( -2.4429 9.4865 -4.5498 \ 
Mij = dB 2 /dxi dB 2 /dx 2 dB 2 /dx 3 = 4.4430 1.4220 2.7926 , (A3) 
V dBs/dx! dB 3 /dx 2 dB 3 /dxs J \ -6.0043 0.6362 0.8396 / 

assuming a length scale of 1 and a unit of Gauss. Note that the local electric current (J) and the Lorentz force (F) 
is completely determined by Mij as well. The trace of Mij is just V • B and should vanish. However, because of 
the linearization (when there might be sub-grid structures) and the computational errors, the zero divergence is not 
strictly satisfied. We estimate the relative error in calculating Mij to be |V • B| / |V x B| = 3.2% (Xiao et al. 2006). 

The behavior of B near the singularity is represented by the three eigenvectors vi , v 2 , and V3 of M^ . We find the 
eigenvectors and their corresponding eigenvalues Ai, A2, and A3 to be: 

Ai = 5.4421, vi = (0.6118, 0.1347, -0.7795) T , 

A 2 = 4.4599, v 2 = (-0.5130, 0.0380, 0.8575) T , (A4) 
A 3 = -10.0833, v 3 = (0.7870, -0.4148, 0.4568) T . 

In this case, all three eigenvalues are real with one negative and two positives. The eigenvector V3 with the sole 
negative eigenvalue A3 determines the initial direction of the "spine" field line. The other two eigenvectors vi and v 2 
define the "fan" plane; whereas the linear combination of them gives the initial directions of the "fan" field lines. The 
fan field lines define the separatrix (fan) surface, which separates different domains of magnetic flux. 

In practice, field lines traced slightly away from the null in the fan plane tend to flow along the separatrix surface. It 
is interesting that vi • v 2 = —0.9771, i.e. they are almost 170° with respect to each other. As a result, the traced field 
lines rapidly converge into two groups, one connecting to Nl, the other N2 (Figure 5(b)), forming a cusp structure that 
looks almost two-dimensional. Further analysis classifies this null point as a positive (fan field lines going outward) 
non-potential null, with current components parallel to the spine and perpendicular to it (Parnell et al. 1996). 

We note that in some frames in the time series, multiple null points appear in adjacent computational cells near 
the modeled loop apexes. Both positive and negative nulls exist in the sample, although the morphology of the closed 
loops remains similar. Such behavior may be related to the uncertainties in modeling and data. More work is needed 
to evaluate the effect of errors on the inferred topology. 



REFERENCES 



Aulanier, G., DeLuca, E. E., Antiochos, S. K., McMullen, R. A., 

& Golub, L. 2000, Astrophys. J., 540, 1126 
Borrero, J. H., Tomczyk, S., Kubo, M., Socas-Navarro, H., Schou, 

J., Couvidat, S., & Bogart, R. 2011, Sol. Phys., 273, 267 
Demoulin, P., Henoux, J. C., Priest, E. R., & Mandrini, C. H. 

1996, Astron. Astrophys., 308, 643 
DeRosa, M. L., et al. 2009, Astrophys. J., 696, 1780 
Filippov, B., Koutchmy, S., & Tavabi, E. 2012, Sol. Phys., in press 
Fletcher, L., Metcalf, T. R., Alexander, D., Brown, D. S., & 

Ryder, L. A. 2001, Astrophys. J., 554, 451 
Forbes, T. G. 2000, J. Geophys. Res., 105, 23153 
Gary, G. A. 2001, Sol. Phys., 203, 71 
Golub, L., et al. 2007, Sol. Phys., 243, 63 

Gopalswamy, N., Makela, P., Xie, H., Akiyama, S., Sz Yashiro, S. 

2009, J. Geophys. Res., 114, 

Greene, J. M. 1992, J. Comput. Phys., 98, 194 

Haynes, A. L., & Parnell, C. E. 2007, Phys. Plasmas, 14, 082107 

— . 2010, Phys. Plasmas, 17, 092903 

Hoeksema, J. T., et al. 2012, Sol. Phys., in preparation 

Howard, R. A., et al. 2008, Space Science Reviews, 136, 67 

Hudson, H. S. 2011, Space Sci. Rev., 158, 5 

Jing, J., Chen, P. F., Wiegelmann, T., Xu, Y., Park, S.-H., & 

Wang, H. 2009, Astrophys. J., 696, 84 
Leka, K. D., Barnes, G., Crouch, A. D., Metcalf, T. R., Gary, 

G. A., Jing, J., & Liu, Y. 2009, Sol. Phys., 260, 83 
Lemen, J. R., et al. 2012, Sol. Phys., 275, 17 
Lin, R. P., et al. 2002, Sol. Phys, 210, 3 
Liu, Y. 2007, Astrophys. J., 654, L171 

Longcope, D. W. 2005, Living Reviews in Solar Physics, 2, 7 
Mandrini, C. H., Demoulin, P., Schmieder, B., Deluca, E. E., 

Pariat, E., & Uddin, W. 2006, Sol. Phys., 238, 293 
Masson, S., Pariat, E., Aulanier, G., & Schrijver, C. J. 2009, 

Astrophys. J., 700, 559 
Metcalf, T. R. 1994, Sol. Phys., 155, 235 

Moore, R. L., Cirtain, J. W., Sterling, A. C, & Falconer, D. A. 

2010, Astrophys. J., 720, 757 

Moreno-Insertis, F., Galsgaard, K., & Ugarte-Urra, I. 2008, 

Astrophys. J., 673, L211 
Pariat, E., Antiochos, S. K., & DeVore, C. R. 2009, Astrophys. J., 

691, 61 



Parnell, C. E., Smith, J. M., Neukirch, T., & Priest, E. R. 1996, 

Phys. Plasmas, 3, 759 
Priest, E., & Forbes, T., eds. 2000, Magnetic reconnection: MHD 

theory and applications 
Reardon, K. P., Wang, Y.-M., Muglach, K., & Warren, H. P. 

2011, Astrophys. J., 742, 119 
Regnier, S., & Canfield, R. C. 2006, Astron. Astrophys., 451, 319 
Sakurai, T. 1989, Space Sci. Rev., 51, 11 
Schou, J., et al. 2012, Sol. Phys., 275, 229 
Schrijver, C. J. 2009, Advances in Space Research, 43, 739 
Schrijver, C. J., Aulanier, G., Title, A. M., Pariat, E., &; 

Delannee, C. 2011, Astrophys. J., 738, 167 
Shibata, K., Masuda, S., Shimojo, M., Hara, H., Yokoyama, T., 

Tsuneta, S., Kosugi, T., &; Ogawara, Y. 1995, Astrophys. J., 

451, L83 

Shibata, K., Nitta, N., Strong, K. T., Matsumoto, R., Yokoyama, 
T., Hirayama, T., Hudson, H., & Ogawara, Y. 1994, Astrophys. 
J., 431, L51 

Shibata, K., Shimojo, M., Yokoyama, T., & Ohyama, M. 1997, in 
Astronomical Society of the Pacific Conference Series, Vol. Ill, 
Magnetic Reconnection in the Solar Atmosphere, ed. R. D. 
Bentley & J. T. Mariska, 29 

Su, J. T., Shen, Y. D., & Liu, Y. 2012, Astrophys. J., 754, 43 

Sun, X., Hoeksema, J. T., Liu, Y., Wiegelmann, T., Hayashi, K., 
Chen, Q., & Thalmann, J. 2012, Astrophys. J., 748, 77 

Thalmann, J. K., & Wiegelmann, T. 2008, Astron. Astrophys., 
484, 495 

Thompson, B., Demoulin, P., Mandrini, C, Mays, M., Ofman, L., 

Van Driel-Gesztelyi, L., & Viall, N. 2011, in AAS/Solar Physics 

Division Abstracts #42, 2117 
Tian, H., Mcintosh, S. W., Xia, L., He, J., & Wang, X. 2012, 

Astrophys. J., 748, 106 
Titov, V. S., Hornig, G., & Demoulin, P. 2002, J. Geophys. Res., 

107, 1164 

Torok, T., Panasenco, O., Titov, V. S., Mikic, Z., Reeves, K. K., 
Velli, M., Linker, J. A., & De Toma, G. 2011, Astrophys. J., 
739, L63 

Tsuneta, S. 1996, Astrophys. J., 456, 840 
Tsuneta, S., et al. 2008, Sol. Phys., 249, 167 



10 



Sun et al. 



Wang, Y., Chen, C, Gui, B., Shen, C, Ye, P., & Wang, S. 2011, Wiegelmann, T., Inhester, B., & Sakurai, T. 2006, Sol. Phys.. 

J. Geophys. Res., 116, 4104 233, 215 

Wiegelmann, T. 2004, Sol. Phys., 219, 87 Xiao, C. J., et al. 2006, Nat. Phys., 2, 478 



